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Abstract 

The structure of the nonperturbative vector meson vertex function complicates the 
proof of the factorisation theorem for the reaction 7*p — > Vp. It leads to additional 
contributions but, in a simple model for the vertex function, gauge invariance 
ensures that they cancel and factorisation is preserved. 



1. Introduction 

With the measurement of exclusive vector meson electroproduction at high photon en- 
ergy and virtuality at HERA [T[| a new possibility for investigating the interplay of hard 
and soft physics is emerging. The original crude calculations of the process ]2| , based on 
the diagrams of figure [I], have since been refined by several authors || . These diagrams 
contain a hard part, calculated with perturbation theory, but there is also a nonpertur- 
bative part. Recently, the mechanism of factorisation of the calculable hard amplitude 
and the nonperturbative component, given by the non-diagonal parton densities and the 
meson wave function, has been discussed in a rather general framework Q. 

It is the purpose of the present paper to discuss the problems associated with the 
structure of the vertex function at the top right-hand corner of the diagrams. This 
couples the vector meson to the internal quark loop. It is nonperturbative and so our 
knowledge of it is far from complete, though its analytic structure is known ||. If we 
neglect spin, it is a function V(u,v) of the squared 4-momenta on the two quark legs, 
u = k 2 and v = (q' — k) 2 . In the case of a nonrelativistic system, such as the T, it is 
a good approximation to take V to be strongly peaked at values of u and v such that 
the quarks are very close to their mass shells. But it is doubtful that this is a good 
approximation for the J/ip, and for the p surely it is not. In this case, V(u, v) has branch 
points in each of its variables, u and v, associated with "normal thresholds" and possibly 
also "anomalous thresholds" f|. We show that even in leading power this results in a 
breakdown of the desired factorisation of the sum of the diagrams of figure |l|. 

However, this does not necessarily imply that the factorisation theorem is invalid. 
The set of diagrams of figure [l] is not QCD gauge invariant. To achieve gauge invariance, 
they must be supplemented by other graphs in which either or both of the gluons couples 
directly into the nonperturbative vector meson vertex function. The question we discuss 
is whether these additional diagrams restore the factorisation theorem. 

We are not able to give a definitive answer to this question, because of the need 
to introduce two further nonperturbative vertex functions, and we know as little about 
these as we do about V. However, we report a calculation based on an explicit simple 
model in which we find that the factorisation theorem is indeed restored. This encourages 
the belief that its validity may be general. 

Our model is a simple one in which all three vertex functions have the analytic struc- 
ture that is expected on general grounds J5J and in which they are related in such a way 
that the complete leading-power amplitude for the exclusive meson production is gauge 
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invariant. We then find that the amplitude contains leading contributions from diagrams 
that cannot be separated into a hard production process and soft wave function correc- 
tions. However, gauge invariance leads to a cancellation of the unwanted contributions 
and leaves the factorisation theorem intact. The final formula can be obtained from the 
diagrams of figure ^ by a redefinition of the rules for calculating them: basically, this 
amounts to ignoring the branch points in the vertex function V. 

In section 2, we review some aspects of the factorisation theorem in the case where 
the structure of the nonperturbative vertex function V may be neglected. In section 3 
we explain the complications that arise when its structure is taken into account, and 
introduce a simple model for this structure. Its consequences are analysed in section 4. 




a) b) c) 

Fig.|] The leading amplitude for a structureless meson vertex. 

2. Structureless vertex function 

According to usual ideas, unless || the vector meson is at least as heavy as the T, its 
complete vertex function V(u, v) is nonperturbative for small values of u and v. Its high-w 
or high-f tail may be calculated from perturbative QCD, but we do not consider this tail 
because it contributes a correction to the amplitude that is a nonleading power of the 
perturbative coupling as{Q 2 )- To begin with, we take the nonperturbative vertex to be 
structureless. For simplicity, we take the meson to be spinless, and also pretend that the 
photon and the quarks have no spin. The coupling of the gluons to the scalar quarks is 
given by — ig (r^+r' ), where r and r' are the momenta of the directed quark lines, and the 
coupling of the photon is ie, where e has the dimensions of mass. We use light-cone co- 
ordinates and work in a frame where q, q' and P, P' have zero transverse components and 
large '+' and ' — ' components respectively. We concentrate on the forward production, 
so that the transverse component of the momentum transfer A = q' — q = £ — £' = P — P' 
vanishes, = 0. 
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For scalar quarks, there is a seagull two-gluon vertex. In order to obtain a gauge- 
invariant set of diagrams, we must supplement those of figure [I] with diagrams that 
involve the seagull. However, we shall calculate in Feynman gauge and suppose that 
W 2 ^> Q 2 , where W is the 7*p energy, W 2 = (q + P) 2 . In this case the diagrams that 
involve the seagull contribute to the amplitude a nonleading power of W: the diagrams 
of figure [I] are the only ones that contribute in leading power when V is constant. We 
study the case of large Q 2 , Q 2 3> m v , so that we neglect all masses. Then for forward 
production A = XBj P. 

The lower bubble in the diagrams of figure |l] in principle has a complicated structure. 
However, we shall assume that the main contribution comes from values of its subenergy 
cr that are not too large, a = (P — t) 2 <C W 2 . Then £ + <C q + and the dependence of 
the lower bubble on a may be approximated by 6(a). Further, when the upper parts 
of the diagrams are added together the main contribution arises from small values of 
£-, L < P_, so that also £'_ <C P_ and I 2 ~ i' 2 ~ £\. This may be seen most simply 
by calculating the imaginary part of the amplitude, where the left-most of the two quarks 
to which £' is attached is on shell. Hence the important part of the lower bubble effectively 
has the structure 

F^(£, £', P) « 5(P„£+) F{£\) P»P U , (1) 

which is defined to include both gluon propagators and all colour factors. A similar 
expression was found by Cheng and Wu [0] in a tree model for the lower bubble, though 
we do not need to restrict ourselves to such a simple model. We assume that F restricts 
the gluon momentum to be soft, £\ <CQ 2 . In the high energy limit it suffices to calculate 

r d A £ r d 4 £ 

M = / -——T^Fuv w / ^-rT ++ F__ , (2) 
J (2vr) 4 J 4(2vr) 4 ++ K J 

where 

T n» = T i*> ^ £f i9 ) = rr + rr + Tp (3) 

is the sum of the upper parts of the diagrams in figure [3]. 

The lower amplitude F^ in the diagrams of figure 1 is symmetric with respect to the 
two gluon lines. This symmetry of the lower amplitude allows us to replace the properly- 
symmetrised upper amplitude Tj^ a with the unsymmetrised amplitude T^ v correspond- 
ing to the sum of the diagrams in figure 1. The symmetrised amplitude is 

TgL(£, £', q) = \[T^{£, £', q) + V"*(-t, -£, q)] . (4) 
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The two exchanged gluons together must form a colour singlet and so, at least for our 
calculation that takes account only of the lowest order in as{Q 2 ), the symmetrised am- 
plitude T^ m satisfies the same Ward identity as for two photons: 

T^ m (£,£',q)£/ u = 0. (5) 

Writing this equation in light-cone components and setting £± = £' ± , we see that for the 
small values of £'_, £ + and £' + that we need, 

T symi++ ~ £]_ (6) 

for £\ — > 0. Here we have used the fact that the tensor T^ u m , which is built from £', £ 
and q, has no large ' — ' components. The £_ integration makes this equation hold also 
for the original, unsymmetrised amplitude: 

f-2 



d£_T ++ ~£ 2 ± . (7) 

This is the crucial feature of the two-gluon amplitude that will simplify the calculation 
and lead to the factorising result of the next sections. 

It is convenient to begin with the contribution from diagram a) of figure [l] to the £_ 
integral of T ++ , which is required in (0): 

d A k z(l-z) V(k 2 ,(q' -k) 2 ) 



J d£_T a , ++ = -4eg 2 q + J 



(2tt) 3 N 2 + (k ± + £ ± ) 2 k 2 (q'-k) 2 ' v/ 

Here N 2 = z(l — z)Q 2 , z = k + /q + and the condition £ + = 0, enforced by the 5-function 
in (0), has been anticipated. 

Now / d£-Tb t++ and / d£-T c<++ each carry no £± dependence. So to ensure the validity 
of (0) the sum of the three diagrams must be 

2 f d 4 k . . . r V(k 2 , (q' - A;) 2 ) 



where 

N 1 1 



(10) 



{N 2 + k\) 2 N 2 + (k ± + £ ± ) 
We have used the softness of the wave function, which results in the dominant contribu- 
tion to the integral arising from values of fc} < Q 2 , and the rotational symmetry, which 
makes k± ■ £± integrate to 0. Note the 1/Q A behaviour obtained after a cancellation of 
1/Q 2 contributions from the individual diagrams. This cancellation, which is closely re- 
lated to the well-known effect of colour transparency ||, has been discussed in in the 
framework of vector meson electroproduction. 
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Introduce the light-cone wave function of the meson 

The final result following from (0) and (|9]) is a convolution of the production amplitude 
of two on-shell quarks and the light-cone wave function: 

M = % eg 2 W 2 ( / ^^1)) JdzJ e^^=^, kl) . (12) 

This corresponds to the 0(£ 2 L ) term in the Taylor expansion of the contribution (f|) from 
figure [pa) . 

In a more general analysis the above £± integral of F has to be replaced by an expres- 
sion proportional to the non-diagonal gluon density which introduces an W 2 dependence 
and an additional Q 2 dependence. 

3. Modelling the meson wave function 

In the previous section we have supposed that the vertex function V has no structure. We 
used the gauge invariance of the amplitude to argue that there is a cancellation among the 
three diagrams of figure |I|a). Of course, we may instead obtain the same result by explicit 
calculation of each of the three diagrams. This may be done most simply by calculating 
their imaginary parts and using the known fact that the complete amplitude must be 
pure-imaginary when its energy dependence is (W 2 ) 1 ' . Alternatively, the contributions 
to the amplitude itself may be calculated by doing the fc_ integration, completing the 
integration contour with an infinite semicircle and taking appropriate pole residues. 

Consider figure |I]b) for example. In either method of calculation, if V has no structure 
the upper quark line gets put on shell. However, if we now take account of the known 
structure of V, there is an additional contribution to the imaginary part corresponding 
to cutting the graph through the vertex function — the k_ integration would need to 
take account also of branch points of V, not just the poles of the propagators. The gauge- 
invariance argument breaks down because the set of diagrams by itself is no longer gauge 
invariant: one must add to it diagrams where either or both of the gluons couples directly 
into the vertex function. 

In order to study this, we use the simplest model for the quark-quark-meson vertex 
function V(u,v) that incorporates at least part of its known branch-point structure. 
It is the purely nonperturbative vertex function that we are modelling: it goes to zero 
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a) 




b) 




Fig.|2]a) Model for the nonperturbative vertex function V. b) Diagram with a gluon 
coupling into the vertex. 

suitably rapidly when either of the squared 4-momenta u or v of the quarks becomes 
large. We do not consider its perturbative tail, which would be obtained by exchanging a 
perturbative gluon between the quarks. It is a familiar notion [0] that the correct analytic 
properties of nonperturbative amplitudes are those corresponding to Feynman graphs, 
even though the numerical values of such graphs have no physical significance. In order 
to model the vertex function V, therefore, we use the simple Feynman graph of figure 0a), 
where the line that joins the quarks is a scalar (like the quarks themselves in our simple 
calculations) which couples to them with strength A' and where the right-hand internal 
vertex is taken to be a constant A. This model has branch points in each of the variables 
u and v, and it has the appropriate softness. 




Fig.|3] Diagram for meson production with the vertex modelled by scalar particle ex- 
change. 

When the vertex function of figure [2]a) is used in figure [T|a), we obtain figure [5| The 
expression for the upper part of the diagram is (|S|) with 

d A k' i\X' 2 



V{k\ (q' - k) ) 



(13) 



(2tt) 4 k' 2 (q' - k') 2 (k - k') 2 ' 
The diagram of figure ^| by itself gives no consistent description of meson production since 
it lacks gauge invariance. This problem is not cured by just adding the two diagrams |T]b) 
and c) with the blob replaced by the vertex of figure 0a). It is necessary also to include 
diagrams where a gluon is coupled into the vertex (see, e.g., figure 0b)). Furthermore, 
diagrams where both gluons couple into the vertex have to be included. 
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As will be demonstrated in the next section, the complete result can nevertheless be 
extracted directly from @ and (|T3|). 



4. Gauge invariance and factorisation 

To obtain gauge invariance within the model of the vertex illustrated in figure 0a), 
diagrams with the two gluons attached to the quark loop in all possible ways have to be 
added to that of figure [| In addition to figure ^| there are nine such diagrams. They are 
shown in figure 




Fig.|4] The remaining diagrams contributing to meson production within the above simple 
model for the meson wave function. 

The same gauge invariance arguments that lead to (|7D apply to the sum of all the 
diagrams in figures [|and|]. Therefore, the complete result for T ++ , which is now defined 
by the sum of the upper parts of all these diagrams, can be obtained by extracting the l\ 
term at leading order in W 2 and Q 2 . Such a term, with a power behaviour ~ £ 2 L W 2 /Q A , 
is obtained from the diagram in figure |3] (see (^) and (|TB|)) by expanding around £±_ = 0. 
We now argue that none of the other diagrams gives rise to such a leading-order l\ 
contribution. 
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It is simpler to evaluate the contributions to the amplitude itself, rather than its 
imaginary part. We continue the evaluation of the diagram in figure |3] by performing the 
fc_ and k'_ integrations in (H) and fll3p . Consider the region where k + > k' + and pick up 
the poles in such a way that the propagators with momenta k! and q' — k go on-shell. 
The resulting expression reads 



-ieg 2 \\' 2 q + f 2 z(l - z) 

' Tfigurc3 ' ++ " (2vr)5 J dzdk± [N 2 + (k ± + £^]kl X (14) 



dz'd 2 k' 



1 



[0 - z'){k'l/z' + fci/(l - z)) + (Aj ± - A;l) 2 ]fc 
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where z' = k' + /q + and the integration is restricted to z > z' . A corresponding expression 
can be obtained for the region z < z' . 

The above manipulations can also be applied to all the diagrams of figure One 
arrives at the following results. The upper parts of diagrams a) - d) have no £± de- 
pendence at all since the £- integration puts the quark propagator connecting the two 
gluon vertices on-shell. Diagrams h) and i) vanish after the £- integration because the 
integrand has two poles, both of which lie on the same side of the real axis. To anal- 
yse the remaining diagrams e) - g) it is convenient to route the momentum £ as far as 
possible to the left. After all the '— ' integrations are performed it becomes obvious that 
each diagram contributes only at order W 2 /Q A . To obtain an £\ term one of the off-shell 
propagators left of the gluon vertices has to be expanded around £± = 0. This brings the 
power down to W 2 /Q 6 . Therefore, no leading-order ^-contribution is obtained. 

The above discussion shows that the complete answer is given by the £\ term from 
the Taylor expansion of (^). The amplitude M is precisely the one of ( |T2| ) and QHD, with 
V(k 2 , (q' - k) 2 ) given by (0). 

We have also checked the correctness of this simple factorising result by explicitly 
reducing all diagrams of figure |to a form similar to that given in (|i4|) for the diagram 



of figure H The required cancellations occur on the level of the integrands, before the 
k±, k' ± , z and z' integrations. 



5. Conclusions 



The mechanism of factorisation in exclusive meson production has been analysed in 
the framework of a simple scalar model. In this model the meson is formed by two 
scalar quarks interacting via the exchange of a scalar boson. From a calculation of all 
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contributing diagrams within the restriction of two-gluon exchange the following picture 
emerges. 

The complete result contains leading contributions from diagrams that cannot be 
factorised into quark-pair production and meson formation. Nevertheless, in Feynman 
gauge the answer to the calculation can be anticipated by looking only at one particular 
factorising diagram. The reason for this simplification is gauge invariance. In the dom- 
inant region where the transverse momentum £± of the two t-channel gluons is small, 
gauge invariance requires the complete quark part of the amplitude to be proportional to 
£ 2 ± . The leading i\ dependence comes exclusively from one diagram. Thus, the complete 
answer can be obtained from this particular diagram, which has the property to factorise 
explicitly if the two quark lines are cut. The resulting amplitude can be written in a 
factorised form. 

It was our intention to demonstrate the above mechanism of gauge invariance induced 
cancellations in as simple and as explicit a way as possible. Our analysis supplements the 
otherwise much more general and complete discussion of [|J by making it more explicit 
and by handling the known structure of the nonperturbative meson vertex function. 

Our calculation has many important limitations. It is restricted to the very simplest 
model of the meson vertex function, with only the simplest features of the known branch- 
point structure. However, a generalisation to a more complicated wave function structure 
seems straightforward. We have assumed that the main contribution to the diagrams in 
figures [TJ, [3] and [| comes from values of the subenergy of the lower bubble that are not too 
large. Nevertheless, the integration over this subenergy will provide the energy growth 
of our cross section, corresponding to the rapid rise at small x of the gluon distribution 
observed at HERA. We have not discussed the subtleties involved in factorising the off- 
diagonal gluon density, nor the behaviour in the region where the longitudinal momentum 
fraction of one of the quarks becomes small. 
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